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ABSTRACT 

We put constraints on the main cosmological parameters of different spatially 
flat cosmological models by combining the recent clustering results of XMM- 
Newton soft (0.5-2 keV) X-ray sources, which have a redshift distribution with 
median redshift 2 ~ 1.2, and SNIa data. Using a likelihood procedure we 
find that the model which best reproduces the observational data and that is 
consistent with stellar ages is the concordance ACDM model with: Q rn ~ 0.28, 
w~ — 1, H 0 — 72 km s -1 Mpc -1 , t Q — 13.5 Gyr and has an X-ray AGN 
clustering evolution which is constant in physical coordinates. For a different 
clustering evolution model (constant in comoving coordinates) we find another 
viable model, although less probable due to the smaller age of the Universe, 
with fi m ~ 0.38, w~ —1.25, H 0 ~ 70 km s -1 Mpc -1 and t 0 ~ 12.9 Gyr. 

Key words: galaxies: clustering- X-ray sources - cosmology:theory - large- 
scale structure of universe 


1 INTRODUCTION 

Recent advances in observational cosmology, based on 
the analysis of a multitude of high quality observa¬ 
tional data (Type la supernovae, cosmic microwave 
background (CMB), large scale structure, age of the 
globular clusters, high redshift galaxies), strongly indi¬ 
cated that we are living in a flat (fi to t = 1) accelerating 
Universe containing a small baryonic component, non- 
baryonic cold dark matter (CDM) to explain the clus¬ 
tering of extragalactic sources and an extra component 
with negative pressure, usually named “dark energy”, 
to explain the present accelerated expansion of the uni¬ 
verse (eg. Riess, et al. 1998; Perlmutter et al. 1999; Ef- 
stathiou et al. 2002; Spergel et al. 2003; Percival et al. 
2002; Tonry et al. 2003; Schuecker et al. 2003; Riess et 
al. 2004; Tegmark et al. 2004) 

The last few years there have been many theoreti¬ 
cal speculations regarding the nature of the exotic “dark 
energy”. Various candidates have been proposed in the 
literature, among which the time varying A-parameter 
(eg. Ozer & Taha 1987), a scalar field having a self¬ 
interaction potential V('I>) with the field energy density 
decreasing with a slower rate than the matter energy 
density (dubbed also “quintessence”, eg. Peebles & Ra- 
tra 2003 and references therein) or an extra “matter” 
component, which is described by an equation of state 
pq = w pQ with w< —1/3 (a redshift dependence of w 


is also possible but present measurments are not pre¬ 
cise enough to allow meaningful constraints; eg. Dicus 
& Repko 2004). A particular case of “dark energy” is the 
traditional A-model which corresponds to w= —1. Note 
that a variety of observations indicate that w< —0.6 for 
a flat geometry (eg. Ettori, Tozzi & Rosati 2002; Tonry 
et al. 2003; Riess et al. 2004; Schuecker 2005). 

In this paper we put constraints on spatially flat 
cosmological models using the recently derived cluster¬ 
ing properties of the XMM -Newton soft (0.5-2 keV) X- 
ray point sources (Basilakos et al. 2005), the SNIa data 
(Tonry et al. 2003) and the age of globular clusters (eg. 
Krauss 2003; Cayrel et al. 2001). Hereafter will use the 
normalized Hubble constant H a = 100 h km s -1 Mpc -1 . 


2 X-RAY AGN CLUSTERING 

In a previous paper (Basilakos et al 2005) we derived 
the angular correlation function of the soft (0.5-2 keV) 
XMM X-ray sources using a shallow (2-10 ksec) wide- 
field survey (~ 2.3 deg. 2 ). A full description of the data 
reduction, source detection and flux estimation are pre¬ 
sented in Georgakakis et al. (2004). Note that the sur¬ 
vey contains 432 point sources within an effective area 
of ~ 2.1 deg 2 (for f x > 2.7x 10 -14 erg cm -2 s -1 ), while 
for fx > 8.8 x 10 -15 erg cm -2 s -1 the effective area of 
the survey is ~ 1.8 deg 2 . In Basilakos et al (2005) we 
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present the details of the correlation function estima¬ 
tion, the various biases that should be taken into ac¬ 
count (the amplification bias and integral constraint), 
the survey luminosity and selection functions as well 
as issues related to possible stellar contamination. In 
particular, the redshift selection function of our X-ray 
sources, derived by using the soft-band luminosity func¬ 
tion of Miyaji, Hasinger & Schmidt (2000), and assum¬ 
ing the realistic luminosity dependent density evolution 
of our sources, predicts a characteristic depth of z — 1.2 
for our sample (for details see Basilakos et al. 2005). 

Our aim here is to compare the theoretical cluster¬ 
ing predictions from different flat cosmological models 
to the actual observed angular clustering of distant X- 
ray AGNs. For the purpose of this study we use Limber’s 
formula which relates the angular, w(9), and the spatial, 
£(r), correlation functions. In the case of a spatially flat 
Universe, Limber’s equation can be written as: 


w{9) = 2 


fo° fo°° sVtogfo z) dxdu 
[J 0 °° x 2 4>(x)dx] 2 


(1) 


where <j>(x) is the selection function (the probability that 
a source at a distance x is detected in the survey) and x 
is the coordinate distance related to the redshift through 


x(z) 



( 2 ) 


with E(z) = [flm (1 + z) 3 + (1 — flm) (1 + z) 13 } 1 ^ 2 and 
j3 = 3(1 + w). The number of objects within a shell 
(z, z + dz) and in a given survey of solid angle ui a is: 


AN 2 
—— = uisX n a 

d z 


<t >{ x ) ( 0 ) E 1 {z) . 


( 3 ) 


where n a is the comoving number density at zero red¬ 
shift. Combining the above system of equations we ob¬ 
tain: 

rr r°° / i j vv 2 

W{e)=2 ^l Gv it) E{z)dz l c(r ’ 2)dM (4) 

where r is the physical separation between two sources, 
having an angular separation, 9 , given by r ~ 
(1 + z ) _1 (u 2 + a : 2 ^ 2 ) 1 ^ 2 (small angle approximation). 
Therefore, in order to estimate the expected w(0) in 
a cosmological model we also need to determine the 
source redshift distribution dN/dz, which as we said 
previously, is estimated by integrating the appropriate 
Miyaji et al. (2000) luminosity function. 


2.1 The Role and Evolution of Galaxy Bias 

It has been claimed that the large scale clustering of 
different mass tracers (galaxies or clusters) is biased 
with respect to the matter distribution (cf. Kaiser 1984; 
Bardeen 1986). It is also an essential ingredient for cold 
dark matter (CDM) models to reproduce the observed 
galaxy distribution (cf. Davis et al. 1985). Within the 
framework of linear biasing (cf. Kaiser 1984; Benson et 
al. 2000 ), the mass-tracer and dark-matter spatial cor¬ 
relations, at some redshift z, are related by: 

£{r,z) = £ DM (r,z)b 2 {z) , (5) 


where b(z) is the bias evolution function. This has been 
shown to be a monotonically increasing function of red¬ 
shift (Mo & White 1996; Matarrese et al. 1997; Basi¬ 
lakos & Plionis 2001 and references therein). Here we use 
the bias evolution model of Basilakos & Plionis (2001, 
2003), which is based on linear perturbation theory and 
the Friedmann-Lemaitre solutions of the cosmological 
held equations. We remind the reader that for the case 
of a spatially hat cosmological model our general bias 
evolution can be written as: 

b(z) = AE(z)+CE(z) j I0^!dy+1 . (6) 

Note that our model gives a family of bias curves, due 
to the fact that it has two unknowns (the integration 
constants A,C). In this paper, for simplicity, we hx the 
value of C to ~ 0.004, as was determined in Basilakos 
& Plionis (2003) from the 2dF galaxy correlation func¬ 
tion. We have tested the robustness of our results by 
increasing C by a factor of 10 and 100 to hnd differences 
of only ~ 5% in the htted values of and b Q . This 
behaviour can be explained from the fact that the dom¬ 
inant term in the right hand side of eq. ( 6 ) is the hrst 
term [cx (1 + z) 3 ^ 2 ] while the second term has a slower 
dependence on redshift [oc (1 + z)\. 


2.2 Clustering Evolution 


The redshift evolution of the spatial mass correlation 
function, (oM(r, z), can be written as the Fourier trans¬ 
form of the spatial power spectrum P(k). Using also eq. 
(5) we have: 


€(r,z) = 


(1 + z) _(3+E V(z) 
2-7T 2 


' fc2p sin(fer) dfc 
kr 


where k is the comoving wavenumber. Note that the pa¬ 
rameter e parametrizes the type of clustering evolution 
(eg. de Zotti et al. 1990). If e = 7 —3 (with 7 the slope of 
the spatial correlation function; 7 = 1 . 8 ) the clustering 
is constant in comoving coordinates while if e = —3 the 
clustering is constant in physical coordinates. 

The power spectrum of our CDM models is given by 
P(k) oc k n T 2 {k) with scale-invariant (n = 1 ) primeval 
inflationary fluctuations and T(k) the CDM transfer 
function. In particular, we use the transfer function pa¬ 
rameterization as in Bardeen et al. (1986), with the cor¬ 
rections given approximately by Sugiyama (1995) while 
the normalization of the power spectrum is given by: 
erg ~ 0.511m 7 with 7 ~ 0.21 — 0.22w + 0.33fl m (Wang & 
Steinhardt 1998). We caution that this fit, based on the 
rich cluster abundances, has been derived for w> —1. In 
this work we assume that the fit is valid also for w< — 1 . 
Note that we also use the non-linear corrections intro¬ 
duced by Peacock & Dodds (1994). 


3 COSMOLOGICAL CONSTRAINTS 

3.1 X-ray AGN Clustering likelihood 

It has been shown that the application of the correlation 
function analysis on samples of high redshift galaxies 
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can be used as a useful tool for cosmological studies (eg. 
Matsubara 2004). In our case to constrain the cosmolog¬ 
ical parameters we utilize a standard x ' 2 likelihood pro¬ 
cedure to compare the measured XMM soft source angu¬ 
lar correlation function (Basilakos et al. 2005) with the 
prediction of different spatially flat cosmological mod¬ 
els. In particular, we define the likelihood estimator* as: 
£ agn (c) oc exp[— Xagn( c )/2] with: 

2 

'tUobs(^i) / 0 \ 

- • 

<7i 

where c is a vector containing the cosmological param¬ 
eters that we want to fit and at the observed angular 
correlation function uncertainty. Here we work within 
the framework of a flat (fi to t = 1 ) cosmology with pri¬ 
mordial adiabatic fluctuations and baryonic density of 
fib h 2 ~ 0.022 (eg. Kirkman et al. 2003). In this case the 
corresponding vector is c = (fi m ,w, h, 6 0 ). We sample 
the various parameters as follows: the matter density 
fi m G [0.01, 1] in steps of 0.01; the equation of state pa¬ 
rameter wG [—3, —0.35] in steps of 0.05, the dimension¬ 
less Hubble constant h G [0.5, 0.9] in steps of 0.02 and 
the X-ray sources bias at the present time b Q G [0.5, 4] 
in steps of 0.05. Note that in order to investigate pos¬ 
sible equations of state, we have allowed the parameter 
w to take values below -1. Such models correspond to 
the so called phantom cosmologies (eg. Caldwell 2002; 
Corasaniti et al. 2004). 

The resulting best fit parameters, for the two clus¬ 
tering evolution models, are presented in the first two 
rows of Table 1. It is important to note that our es¬ 
timates for the Hubble parameter h are in very good 
agreement with those derived (h = 0.72 ± 0.07) by the 
HST key project (Freeman et al. 2001). In Fig.l we 
present the la, 2a and 3a confidence levels in the (w, h ) 
and (w, b 0 ) planes by marginalizing the former over fi m 
and b 0 and the latter over fi m and h. It is evident that 
w is degenerate with respect to both h and the bias at 
the present time. 

Figure 2 also shows the la, 2a and 3a confi¬ 
dence levels (continuous lines) in the (fi m ,w) plane by 
marginalizing over the Hubble constant and the bias 
factor at the present time f The equation of state pa¬ 
rameter likelihood is not constrained by this analysis 
and all the values in the interval —3 < w < —0.35 are 
acceptable within the la uncertainty. Therefore, in or¬ 
der to put further constraints on w we additionally use 
the SNIa data (Tonry et al. 2003) as well as the so called 
age limit, given by the age of the oldest globular clusters 
in our Galaxy (to > 12.7 Gyr; Krauss 2003; Cayrel et 
al. 2001 and references therein). 


Xagn(c) — 


3.2 The likelihood from the SNIa 

We use the sample of 172 supernovae of Tonry et al. 
(2003) in order to constrain and the equation of 


* Likelihoods are normalized to their maximum values, 
t Hereafter, when we marginalize over the Hubble constant 
we will use h = 0.72 for e = —1.2 and h = 0.7 for e = —3. 




w w 

Figure 1. Likelihood contours in the (w, h) plane (left panel) 
and the (w, b 0 ) plane (right panel) for e = —1.2 (a similar 
degeneracy is true also for the e = — 3 clustering evolution 
model). The contours are plotted where —21n£/£ ma x is equal 
to 2.30, 6.16 and 11.83, corresponding to lcr, 2cr and 3cr con¬ 
fidence level. 


state in the framework of a flat geometry (fitot = 1). In 
this case, the corresponding vector c is: c = (fl m , w) and 
the likelihood function can be written as: £ SNIa (c) oc 
ex P[— xi N ia( c )/2] With: 


XSNIa( c ) 


172 

E 


i=1 



(9) 


where Dl(«) is the dimensionless luminosity distance, 
Dh(z) = H 0 (l + z)x(z) and z% is the observed redshift. 
The green lines in Fig. 2 represents the la, 2a,and 3a, 
confidence levels in the (fi m , w) plane. We find that the 
best fit solution is fi m = 0.30 ± 0.04 for w> —1, in 
complete agreement with previous SNIa studies (Tonry 
et al. 2003; Riess et al. 2004). 


3.3 The joined likelihoods 

In order to combine the X-ray AGN clustering proper¬ 
ties with the SNIa data we perform a joined likelihood 
analysis and marginalizing the X-ray clustering results 
over h and b a , which are not constrained by the value of 
w (see Fig.l), we obtain: £ joint (fi m , w) = £ AGN x £ SNIa . 
Taking also into account the age limit (the yellow area 
in Fig. 2) the likelihood for the e = —1.2 cluster¬ 
ing evolution model peaks at fi m = 0.28 ± 0.02 with 
w= —1.05 A q' 2 (corresponding to to = 13.5 Gyr) which 
is in excellent agreement with the WMAP results of 
Spergel et al. (2003) and the REFLEX X-ray clusters + 
SNIa results of Schuecker et al. (2003). For the e = —3 
clustering evolution model we obtain fi m = 0.38 ± 0.03 
with w= —1.25jfQ 25 (which corresponds to to = 12.9 
Gyr). The latter model appears to be marginally ruled 
out by the stellar ages. Note that the normalization of 
the power spectrum that corresponds to these models is 
as — 0.98 and 0.90, respectively. 

It is evident that the combined likelihood analysis 
puts strong constraints on the value of w and once in¬ 
cluding stellar ages it appears to favor the standard con¬ 
cordance ACDM (fi ~ 0.3, w~ —1) cosmological model 
as well as a comoving AGN clustering model (e = —1.2). 
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Table 1. Cosmological parameters from the likelihood analysis: The 1 st column indicates the data used (the last two rows 
correspond to the joint likelihood analysis). Errors of the fitted parameters represent lcr uncertainties. Note that for the joined 
analysis the corresponding results are marginalized over the Hubble constant and the bias factor at the present time, for which 
we use the values indicated. 


Data 


w 

h 

bo 

to 

X 2 /dof 

XMM (e = -1.2) 
XMM(e = -3.0) 

0 31+ 016 
u '° -0.08 
n OO+0.02 
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r» 79+0.02 
u - ‘ z -0.18 

0 70 +0 04 
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2 30 +0 ' 70 
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1 20 +0 ‘ 60 
1-zu —0.30 
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0.84 

XMM(e = —1.2)/SNIa 

0.28 ± 0.02 

— 1 05+ 010 

l.UO — o 20 

0.72 

2.30 

13.5 

0.87 

XMM(e = —3.0)/SNIa 

0.38 ± 0.03 

1 or+0.10 

i.zO—o 25 

0.70 

1.20 

12.9 
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Figure 2. Likelihood contours in the (f2 m ,w) plane. The 
contours correspond to la (2.30), 2 a (6.16) and 3<r (11.83) 
confidence levels, using the two different clustering behaviors 
(left panel for e = —3 and right panel for e = —1.2). Note that 
the black and the green lines correspond to the X-ray AGN 
clustering and SNIa results, respectively while the yellow area 
is ruled out by the stellar ages. 


However, the model with w~ —1.25, — 0.38 and 

e = —3 cannot be ruled out at any significant level. 

Many other recent analyses utilizing different com¬ 
binations of data seem to agree with the former cosmo¬ 
logical model. For example, Tegmark et al. (2004) used 
the WMAP CMB anisotropies in combination with the 
SDSS galaxy power spectrum and found a good ACDM 
fit with = 0.30 ± 0.04 and h = 0.70j(QQ3 (see also 
Spergel et al. 2003; Percival et al. 2003; Schuecker et al. 
2003). Also combining the gas fraction in relaxed X-ray 
luminous clusters with the CMB and SNIa has provided 
stringent constrains with ~ 0.3 and w ~ — 1 (eg. 
Allen et al. 2004; Rapetti, Allen & Weller 2004). 


4 CONCLUSIONS 

We have combined for the first time the clustering prop¬ 
erties of distant X-ray AGNs, identified as soft (0.5-2 
keV) point sources in a shallow ~ 2.3 deg 2 XMM sur¬ 
vey, which have a 2 -distribution that peaks at 2 ~ 1 . 2 , 
with the SNIa data. From the X-ray AGN clustering 
likelihood analysis alone we constrain h ~ 0.72 ± 0.03 
(where the uncertainty is found after marginalizing over 
w and b 0 ). From the joined likelihood analysis and tak¬ 
ing into account stellar ages we constrain the matter 
density and the equation of state parameters. The best 
model appears to be one with S7 m ~ 0.28, w~ —1 and 


a stable in comoving coordinates X-ray AGN clustering 
model. However, the model with ~ 0.38, w~ —1.25 
and constant in physical coordinates (e = —3) X-ray 
AGN clustering, of which the predicted age is marginally 
consistent with stellar ages, cannot be excluded at any 
significant level by our analysis. 
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ABSTRACT 

We put constraints on the main cosmological parameters of different spatially 
flat cosmological models by combining the recent clustering results of XMM- 
Newton soft (0.5-2 keV) X-ray sources, which have a redshift distribution with 
median redshift z ~ 1.2, and SNIa data. Using a likelihood procedure we 
find that the model which best reproduces the observational data and that is 
consistent with stellar ages is the concordance ACDM model with: f2 m ~ 0.28, 
w~ —1, H 0 ~ 72 km s -1 Mpc -1 , t, 0 — 13.5 Gyr and has an X-ray AGN 
clustering evolution which is constant in physical coordinates. For a different 
clustering evolution model (constant in comoving coordinates) we find another 
viable model, although less probable due to the smaller age of the Universe, 
with fl m ~ 0.38, w~ —1.25, H a ~ 70 km s -1 Mpc -1 and t 0 — 12.9 Gyr. 

Key words: galaxies: clustering- X-ray sources - cosmology:theory - large- 
scale structure of universe 


1 INTRODUCTION 

Recent advances in observational cosmology, based on 
the analysis of a multitude of high quality observa¬ 
tional data (Type la supernovae, cosmic microwave 
background (CMB), large scale structure, age of the 
globular clusters, high redshift galaxies), strongly indi¬ 
cated that we are living in a flat (fi to t = 1) accelerating 
Universe containing a small barvonic component, non- 
baryonic cold dark matter (CDM) to explain the clus¬ 
tering of extragalactic sources and an extra component 
with negative pressure, usually named “dark energy”, 
to explain the present accelerated expansion of the uni¬ 
verse (eg. Riess, et al. 1998; Perlmutter et al. 1999; Ef- 
statliiou et al. 2002; Spergel et al. 2003; Percival et al. 
2002; Tonry et al. 2003; Schuecker et al. 2003; Riess et 
al. 2004; Tegmark et al. 2004) 

The last few years there have been many theoreti¬ 
cal speculations regarding the nature of the exotic “dark 
energy”. Various candidates have been proposed in the 
literature, among which the time varying A-parameter 
(eg. Ozer & Talia 1987), a scalar field having a self- 
interaction potential V(4>) with the field energy density- 
decreasing with a slower rate than the matter energy 
density (dubbed also “quintessence”, eg. Peebles & Ra- 
tra 2003 and references therein) or an extra “matter” 
component, which is described bv an equation of state 
PQ = w pq with w< —1/3 (a redshift dependence of w 


is also possible but present measurments are not pre¬ 
cise enough to allow meaningful constraints; eg. Dicus 
& Repko 2004). A particular case of “dark energy” is the 
traditional A-model which corresponds to w= —1. Note 
that a variety of observations indicate that w< —0.6 for 
a flat geometry (eg. Ettori, Tozzi & Rosati 2002; Tonry 
et al. 2003; Riess et al. 2004; Schuecker 2005). 

In this paper we put constraints on spatially flat 
cosmological models using the recently derived cluster¬ 
ing properties of the XMM -Newton soft (0.5-2 ke\' T ) X- 
ray point sources (Basilakos et al. 2005), the SNIa data 
(Tonry et al. 2003) and the age of globular clusters (eg. 
Krauss 2003; Cayrel et al. 2001). Hereafter will use the 
normalized Hubble constant H 0 = 100 hkms -1 Mpc -1 . 


2 X-RAY AGN CLUSTERING 

In a previous paper (Basilakos et al 2005) we derived 
the angular correlation function of the soft (0.5-2 keV) 
XMM X-ray sources using a shallow (2-10ksec) wide- 
field survey (~ 2.3 deg. 2 ). A full description of the data 
reduction, source detection and flux estimation are pre¬ 
sented in Georgakakis et al. (2004). Note that the sur¬ 
vey contains 432 point sources within an effective area 
of ~ 2.1 deg 2 (for f x > 2.7x 10 -14 erg cm -2 s -1 ), while 
for f x > 8.8 x 10 -1 ° erg cm -2 s -1 the effective area of 
the survey is ~ 1.8 deg 2 . In Basilakos et al (2005) we 
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present the details of the correlation function estima¬ 
tion, the various biases that should be taken into ac¬ 
count (the amplification bias and integral constraint), 
the survey luminosity and selection functions as well 
as issues related to possible stellar contamination. In 
particular, the redsliift selection function of our X-ray 
sources, derived by using the soft-band luminosity func¬ 
tion of Miyaji, Hasinger & Schmidt (2000), and assum¬ 
ing the realistic luminosity dependent density evolution 
of our sources, predicts a characteristic depth of z — 1.2 
for our sample (for details see Basilakos et al. 2005). 

Our aim here is to compare the theoretical cluster¬ 
ing predictions from different flat cosmological models 
to the actual observed angular clustering of distant X- 
ray AGNs. For the purpose of this study we use Limber’s 
formula which relates the angular, w(8), and the spatial, 
£(r), correlation functions. In the case of a spatially flat 
Universe, Limber’s equation can be written as: 


w{9) = 2 


L 00 JT x 4 <f> 2 (x.)£,(r, z)dr:dn 
[/ 0 °° x 2 (/){x)dx]‘ 2 


( 1 ) 


where cj)(x) is the selection function (the probability that 
a source at a distance x is detected in the survey) and x 
is the coordinate distance related to the redsliift through 


x(z) 


c f Z dy 
HoJ o E(y) 


( 2 ) 


with E{z) = [fi m (l + z) 3 + (1 — fl m )(l + zf] 1/2 and 
= 3(1 + w). The number of objects within a shell 
(z, z + dz) and in a given survey of solid angle w s is: 


d N 
dz 


= Lu s x 2 n s </>(x) E 1 (z) . 


(3) 


where n s is the comoving number density at zero red- 
shift. Combining the above system of equations we ob¬ 
tain: 

- = (vf) V-)d*£“ { (r,.-)d„ (4) 

where r is the physical separation between two sources, 
having an angular separation, 9, given by r ~ 
(i + sr 1 (m 2 + x 2 e 2 ) 1/ " (small angle approximation). 
Therefore, in order to estimate the expected w(8) in 
a cosmological model we also need to determine the 
source redsliift distribution dN/dz , which as we said 
previously, is estimated by integrating the appropriate 
Miyaji et al. (2000) luminosity function. 


2.1 The Role and Evolution of Galaxy Bias 

It has been claimed that the large scale clustering of 
different mass tracers (galaxies or clusters) is biased 
with respect to the matter distribution (cf. Kaiser 1984; 
Bardeen 1986). It is also an essential ingredient for cold 
dark matter (CDM) models to reproduce the observed 
galaxy distribution (cf. Davis et al. 1985). Within the 
framework of linear biasing (cf. Kaiser 1984; Benson et 
al. 2000), the mass-tracer and dark-matter spatial cor¬ 
relations, at some redsliift z, are related by: 

£0,2) = £ D M(r,z)& 2 (z) , (5) 


where b(z) is the bias evolution function. This has been 
shown to be a monotonically increasing function of red- 
shift (Mo & White 1996; Matarrese et al. 1997; Basi¬ 
lakos & Plionis 2001 and references therein). Here we use 
the bias evolution model of Basilakos & Plionis (2001, 
2003), which is based on linear perturbation theory and 
the Friedmann-Lemaitre solutions of the cosmological 
field equations. We remind the reader that for the case 
of a spatially flat cosmological model our general bias 
evolution can be written as: 

b(z) = AE(z) + CE(z) j ^Jt^dy + l . (6) 

Note that our model gives a family of bias curves, due 
to the fact that it has two unknowns (the integration 
constants A,C). In this paper, for simplicity, we fix the 
value of C to ~ 0.004, as was determined in Basilakos 
& Plionis (2003) from the 2dF galaxy correlation func¬ 
tion. We have tested the robustness of our results by 
increasing C by a factor of 10 and 100 to find differences 
of only ~ 5% in the fitted values of f2 m and b 0 . This 
behaviour can be explained from the fact that the dom¬ 
inant term in the right hand side of eq. (6) is the first 
term [oc (1 + z) 3 ^ 2 ] while the second term has a slower 
dependence on redsliift [oc (1 -f z)]. 


2.2 Clustering Evolution 


The redsliift evolution of the spatial mass correlation 
function, £dm(u z), can be written as the Fourier trans¬ 
form of the spatial power spectrum P(k). Using also eq. 
(5) we have: 


£(r,z) 


(i + zr ( 3 +e, & 2 (z) r 

2^’ Jo 


k 2 P(k) 


., 7 , 


where k is the comoving wavenumber. Note that the pa¬ 
rameter e parametrizes the type of clustering evolution 
(eg. de Zotti et al. 1990). If e = 7 —3 (with 7 the slope of 
the spatial correlation function; 7 = 1.8) the clustering 
is constant in comoving coordinates while if e = —3 the 
clustering is constant in physical coordinates. 

The power spectrum of our CDM models is given by 
P(k) oc k n T 2 (k) with scale-invariant (n = 1) primeval 
inflationary fluctuations and T(k) the CDM transfer 
function. In particular, we use the transfer function pa¬ 
rameterization as in Bardeen et al. (1986), with the cor¬ 
rections given approximately by Sugiyama (1995) while 
the normalization of the power spectrum is given by: 
as — 0.5S2...7 with 7 ~ 0.21 — 0.22w + 0.33Q m (Wang & 
Steinhardt 1998). We caution that this fit, based on the 
rich cluster abundances, has been derived for w> —1. In 
this work we assume that the fit is valid also for w< —1. 
Note that we also use the non-linear corrections intro¬ 
duced by Peacock & Dodds (1994). 


3 COSMOLOGICAL CONSTRAINTS 

3.1 X-ray AGN Clustering likelihood 

It has been shown that the application of the correlation 
function analysis on samples of high redsliift galaxies 
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can be used as a useful tool for cosmological studies (eg. 
Matsubara 2004). In our case to constrain the cosmolog¬ 
ical parameters we utilize a standard \ 2 likelihood pro¬ 
cedure to compare the measured XMM soft source angu¬ 
lar correlation function (Basilakos et al. 2005) with the 
prediction of different spatially flat cosmological mod¬ 
els. In particular, we define the likelihood estimator* as: 
£ agn (c) oc exp[—\' a G n( c )/ 2 ] with: 

Tt r -I 2 

2 , \ _sr^ wth(#*, c) — w obs (0i) 

Xagn( c ) — > - • (°) 

G i 

; = 1 L J 

where c is a vector containing the cosmological param¬ 
eters that we want to fit. and <r, the observed angular 
correlation function uncertainty. Here we work within 
the framework of a flat (fitot = 1 ) cosmology with pri¬ 
mordial adiabatic fluctuations and baryonic density of 
Ob h 2 — 0.022 (eg. Kirkman et al. 2003). In this case the 
corresponding vector is c = (O m ,w ,h,b 0 ). We sample 
the various parameters as follows: the matter density 
O m E [ 0 . 01 , 1 ] in steps of 0 . 01 ; the equation of state pa¬ 
rameter wE [—3, —0.35] in steps of 0.05, the dimension¬ 
less Hubble constant h E [0.5, 0.9] in steps of 0.02 and 
the X-ray sources bias at the present time b a E [0.5, 4] 
in steps of 0.05. Note that in order to investigate pos¬ 
sible equations of state, we have allowed the parameter 
w to take values below -1. Such models correspond to 
the so called phantom cosmologies (eg. Caldwell 2002 ; 
Corasaniti et al. 2004). 

The resulting best fit parameters, for the two clus¬ 
tering evolution models, are presented in the first two 
rows of Table 1. It is important to note that our es¬ 
timates for the Hubble parameter h are in very good 
agreement with those derived (h = 0.72 ± 0.07) by the 
HST key project (Freeman et al. 2001). In Fig.l we 
present the ler, 2<r and 3<r confidence levels in the (w, h) 
and (w, bo) planes by marginalizing the former over Q m 
and bo and the latter over fi m and h. It is evident that 
w is degenerate with respect to both h and the bias at 
the present time. 

Figure 2 also shows the la, 2<r and 3cr confi¬ 
dence levels (continuous lines) in the (fl m ,w) plane by 
marginalizing over the Hubble constant and the bias 
factor at the present time t. The equation of state pa¬ 
rameter likelihood is not constrained by this analysis 
and all the values in the interval —3 < w < —0.35 are 
acceptable within the lcr uncertainty. Therefore, in or¬ 
der to put further constraints on w we additionally use 
the SNIa data (Tonry et al. 2003) as well as the so called 
age limit, given by the age of the oldest globular clusters 
in our Galaxy (to > 12.7 Gyr; Krauss 2003; Cayrel et 
al. 2001 and references therein). 

3.2 The likelihood from the SNIa 

We use the sample of 172 supernovae of Tonry et al. 
(2003) in order to constrain Q m and the equation of 

* Likelihoods are normalized to their maximum values. 

Hereafter, when we marginalize over the Hubble constant 
we will use h = 0.72 for e = —1.2 and h = 0.7 for e = —3. 



Figure 1 . Likelihood contours in the (w ,h) plane (left panel) 
and the (w, bo) plane (right panel) for e = —1.2 (a similar 
degeneracy is true also for the e = — 3 clustering evolution 
model). The contours are plotted where —21n£/£ max is equal 
to 2.30, 6.16 and 11.83, corresponding to lcr, 2<r and 3cr con¬ 
fidence level. 


state in the framework of a flat geometry (Qtot = 1). In 
this case, the corresponding vector c is: c = (Q m , w) and 
the likelihood function can be written as: £ bNIa (c) oc 
ex P[—,\'i N i a ( c )/2] with: 



where £>l(s) is the dimensionless luminosity distance, 
Dl(s) = H 0 (l + z)x(z) and s, is the observed redshift. 
The green lines in Fig. 2 represents the lcr, 2<r,and 3cr, 
confidence levels in the (fi m , w) plane. We find that the 
best fit solution is Q m = 0.30 ± 0.04 for w> —1, in 
complete agreement with previous SNIa studies (Tonry 
et al. 2003; Riess et al. 2004). 


3.3 The joined likelihoods 

In order to combine the X-ray AGN clustering proper¬ 
ties with the SNIa data we perform a joined likelihood 
analysis and marginalizing the X-ray clustering results 
over h and b 0 , which are not constrained by the value of 
w (see Fig.l), we obtain: £ jomt (fi m , w) = £ AGN x £ bNIa . 
Taking also into account the age limit (the yellow area 
in Fig. 2) the likelihood for the e = —1.2 cluster¬ 
ing evolution model peaks at Q m = 0.28 ± 0.02 with 
w= —1.05lg o (corresponding to to = 13.5 Gyr) which 
is in excellent agreement with the WMAP results of 
Spergel et al. (2003) and the REFLEX X-ray clusters + 
SNIa results of Schuecker et al. (2003). For the e = — 3 
clustering evolution model we obtain fl m = 0.38 ± 0.03 
with w= —1.25ig 25 (which corresponds to to = 12.9 
Gyr). The latter model appears to be marginally ruled 
out by the stellar ages. Note that the normalization of 
the power spectrum that corresponds to these models is 
<Ts — 0.98 and 0.90, respectively. 

It is evident that the combined likelihood analysis 
puts strong constraints on the value of w and once in¬ 
cluding stellar ages it appears to favor the standard con¬ 
cordance ACDM (fl ~ 0.3, w~ —1) cosmological model 
as well as a comoving AGN clustering model (e = —1.2). 
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Table 1. Cosmological parameters from the likelihood analysis: The 1 st column indicates the data used (the last two rows 
correspond to the joint likelihood analysis). Errors of the fitted parameters represent lcr uncertainties. Note that for the joined 
analysis the corresponding results are marginalized over the Hubble constant and the bias factor at the present time, for which 
we use the values indicated. 


Data 

Om 

w 

h 

bo 

to 

X 2 /dof 

XMM (e = -1.2) 
XMM(e = -3.0) 

0-3l4o.J 6 8 

°- 38 ±o.?4 

uncons. (w= —1) 
uncons. (w= —1) 

n 70+0.02 
U *' Z -0.18 

0 70+ 0-04 
U, ' U -0.16 

2 30+ 0 - 70 

0.20 

1 20”*”°' 60 
i,ZU — 0.30 

13.0 

12.6 

0.82 

0.84 

XMM(e = —1.2)/SNIa 

0.28 ± 0.02 

— 1 O'r - * -010 

1,uo — 0.20 

0.72 

2.30 

13.5 

0.87 

XMM( £ = —3.0)/SNIa 

0.38 dfc 0.03 

i oc+0.10 

1 - ZO -0.25 

0.70 

1.20 

12.9 

0.85 



0.3 0.4 0.6 0.8 0.3 0.4 0.6 0.8 

n - 

Figure 2. Likelihood contours in the (Q m ,w) plane. The 
contours correspond to lcr (2.30), 2cr (6.16) and 3cr (11.83) 
confidence levels, using the two different clustering behaviors 
(left panel for e = — 3 and right panel for e = —1.2). Note that 
the black and the green lines correspond to the X-ray AGN 
clustering and SNIa results, respectively while the yellow area 
is ruled out by the stellar ages. 

However, the model with w~ —1.25, fl m — 0.38 and 
e = —3 cannot be rilled out at any significant level. 

Many other recent analyses utilizing different com¬ 
binations of data seem to agree with the former cosmo¬ 
logical model. For example, Tegmark et al. (2004) used 
the WMAP CMB anisotropies in combination with the 
SDSS galaxy power spectrum and found a good ACDM 
fit with n m = 0.30 ± 0.04 and h = 0.704^o3 (see also 
Spergel et al. 2003; Percival et al. 2003; Schuecker et al. 
2003). Also combining the gas fraction in relaxed X-ray 
luminous clusters with the CMB and SNIa has provided 
stringent constrains with O m ~ 0.3 and w ~ — 1 (eg. 
Allen et al. 2004; Rapetti, Allen & Weller 2004). 


4 CONCLUSIONS 

We have combined for the first time the clustering prop¬ 
erties of distant X-ray AGNs, identified as soft (0.5-2 
keV) point sources in a shallow ~ 2.3 deg 2 XMM sur¬ 
vey, which have a s-distribution that peaks at 2 ~ 1.2, 
with the SNIa data. From the X-ray AGN clustering 
likelihood analysis alone we constrain h ~ 0.72 ± 0.03 
(where the uncertainty is found after marginalizing over 
w and b 0 ). From the joined likelihood analysis and tak¬ 
ing into account stellar ages we constrain the matter 
density and the equation of state parameters. The best 
model appears to be one with fi m ~ 0.28, w~ —1 and 


a stable in comoving coordinates X-ray AGN clustering 
model. However, the model with ~ 0.38, w~ —1.25 
and constant in physical coordinates (e = —3) X-ray 
AGN clustering, of which the predicted age is marginally 
consistent with stellar ages, cannot be excluded at any 
significant level by our analysis. 
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